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Abstract 



The analysis of complex reaction networks is of great importance in several chemical and biochemical fields (inter- 
stellar chemistry, prebiotic chemistry, reaction mechanism, etc). In this article, we propose to simultaneously refine 
and extend for general chemical reaction systems the formalism initially introduced for the description of metabolic 
networks. The classical approaches through the computation of the right null space leads to the decomposition of 
the network into complex "cycles" of reactions concerned with all metabolites. We show how, departing from the 
left null space computation, the flux analysis can be decoupled into linear fluxes and single loops, allowing a more 
refine qualitative analysis as a function of the antagonisms and connections among these local fluxes. This analysis is 
made possible by the decomposition of the molecules into elementary subunits, called "reactons" and the consequent 
decomposition of the whole network into simple first order unary partial reactions related with simple transfers of 
reactons from one molecule to another. This article explains and justifies the algorithmic steps leading to the total 
decomposition of the reaction network into its constitutive elementary subpart. 
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Introduction 

The dynamical analysis of complex reaction networks such as the ones characterizing interstellar chemistry (Q]), com- 
plex reaction mechanism (2) or metabolic functions (3) should be facilitated by algorithmic means to decouple these 
networks. It is conceivable that the presence of interesting dynamical phenomena like bifurcation or symmetry break- 
ing is mainly due to structural antagonisms between reaction sub-networks provoking "threshold" effects. We propose 
in this article various algorithmic tools to allow such decoupling of complex reaction networks into simpler sub- 
networks. Any of these sub-networks will be restrictively concerned with the successive transformations of one given 
chemical group that will carry the name of "reacton" (in reminiscence of the organic chemistry "synthon"), defined as 
parts of the molecules that are never broken into smaller pieces by any reaction of the network (but internal rearrange- 
ments of the reactons are possible). 

There has been a great amount of literature dedicated to the analysis of these reaction networks into sets of balanced 
reaction fluxes obtained by computing the right null space (or the "row" null space) of the stoichiometric matrix (Q). 
However such analysis, while being performed on the complete molecules, leads to the discovery of fluxes that cross 
the whole set of reactions, making difficult the detection of antagonistic subsets and the automatic anticipation of 
interesting dynamical phenomena. 

The computing of the left null space has been investigated by Famili and Palsson (5). They have shown that 
this leads to "pools of conserved metabolites" through the reaction network. Similar descriptions lead to concepts of 
"conservation analysis" (6) or "metabolic flux analysis" (Q). Beyond what has been proposed in these previous studies, 
this leads to the total decomposition of each molecule into elementary subunits, the reactons. We argue later in the 
article that, such a conservation analysis allows to simplify the global reaction network into smaller sub-networks. 

In order to simplify the reaction network analysis, the first step proposed in this article is to compute the left null 
space (or the "column" null space). The investigation of an optimal basis of this null space (this notion of "optimal 
basis" will be clarified later) reflects the existence of elementary reactons. Once the reactons are identified, the second 
step is to decouple the complete reaction network by simply restricting each sub-network to the single reacton that it 
is concerned with. Provided the reacton basis is correctly selected, it leads to substantially smaller networks. 

Taking for instance the following reaction network: 
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the result of the algorithmic analysis performed on this system should automatically lead to the identification of an 
autocatalytic cycle based one the X%, X3 and Xq compounds, as they are being recycled, of two incoming fluxes of 
X\ and X4, and two outgoing fluxes of X5 and X 2 , as represented in Figure Q] It is easy to detect the presence of two 
linear fluxes and one cycle, the most convenient building blocks to perform the following dynamical analysis (cycles 
and their intrinsic positive feedback are responsible for interesting dynamical effects) and to allow the anticipation of 
any interesting out-of-equilibrium phenomena such as the appearance of bifurcations. 

1 Definition of the network 

The studied chemical network is composed of n different compounds noted for 1 < i < n, and r different 
transformations noted Rj for 1 < j < r. All these reactions are complete chemical transformations, and are thus mass 
balanced between reactant and products. In the example system of Eq.Q] we have n = 6 and r = 3. 
Each reaction can be written in the form: 

n 

Vj€[l,r] Rj : £ v^X % = 0. (4) 

i=i 

vij is the stoichiometric coefficient of the compound i in the transformation j. It is by convention positive for products 
(formed compounds), and negative for reactants (disappearing compounds). 



These stoichiometric coefficients can be gathered in a single n x r matrix v . In this work, all matrices will be 
noted with a *~ . The matrix for the system given in the example is the following 6x3 matrix: 



/ 1 

I 


U 


n \ 
1 


1 

— 1 


n 
U 


1 o 




i 

— i 


n 

u 





-i 








+i 





\o 


+i 


-V 



(5) 



This stoichiometric matrix is actually formed by the juxtaposition of all reaction column vectors: 

V = ("rMT^I ■ ■ ■ |~Tr) . 



(6) 



All column vectors will be noted with a ~~ *, to be distinguished from the row vectors noted with a *~. The purpose 
of this distinction is to easily know the dimension of the vectors, ~a being of dimension n X 1, and "a" of dimension 
1 x r. 

By noting X as the column vector of the JQ, Eq.[4]can be written as a matrix multiplication: 



or, separating each reaction: 
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The traditional analysis of the stoichiometric matrix relies on the calculation of its null spaces. The right null 
space, or the set of combinations of columns of the matrix that gives can be expressed as Q): 
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This conesponds to the set of cycles of reactions, i.e. the combinations of reactions that results in no global change 
inside the system (0). A base of this null space can be computed by the Gauss- Jordan elimination (@), and represented 
by: 
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All the o\ are linearly independent row vectors of the null space, each one being the row k of the matrix The total 
matrix a is thus full rank, and all vectors of the null space can be expressed as a linear combination of the a^. 

On the other hand, the left null space, i.e. the set of combinations of rows of the matrix that gives , can be 
expressed as (@): 
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This left null space of T/ is the right null space of T/ T . It corresponds to the mass balance in compounds (qI). 
Similarly, a base of this null space can be represented as: 



NuU( v) : a = (al\ 
where at indicates the column k of the matrix a . 
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Molecule Decomposition into reactons 



Elementary Reactons 

A reacton is defined as a subpart of a molecule that is never broken into smaller parts by any of the reactions composing 
the network. The chemical reaction network can thus be seen as simple recombinations of reactons. 

A first obvious category of reactons is composed of the atoms, but larger groups are more likely. Typically, in a 
polymerization system, monomers could be such reactons. 

Molecules can be seen as linear combinations of reactons such as: 

a 

Xi = ^2a ijk A k (14) 

k=l 

44> X = *a~A, (15) 

with a being the total number of different reactons chosen for the decomposition. All the are positive or zero 
integers, and represent the number of reacton A k in Xj. The decomposition of Xi is represented by the vector of. The 
reactons can be represented as a column vector A of dimension a. The matrix a of dimension n x a represents the 
decomposition of the whole set of molecules. 

As molecules can always be at least decomposed into unbreakable atoms (no nuclear reaction is obviously consid- 
ered here), there always exists at least one possible combination of reactons. Nevertheless, the computation of the left 
null space should provide us with more useful reactons in order to simplify the whole network into various subsets of 
"independent" reactive pathways. 

Any reaction can be written in terms of reactons by substituting Eq.[l4]in Eq.@} 

n a 
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A T *a T *v' = (18) 



There is a mass balance between each reacton A k : they are never broken into smaller compounds, so that there is 
always the same number of every reacton for both reactants and products in any reaction. We thus have: 



Vj E [1, r] , Vfc € [1, a] , ^Vi,ja ijk = (19) 

i=l 

^ VJfc€[l,a] , a*; T V = 0" (20) 
Thus, we have according to Eq.[[2]and Eq.[T9l 



a k G Null(V) (21) 

All the vectors describing the decomposition of a molecule X{ into reactons A k , as it can be easily obtained by the 
atomic decomposition, are in the left null space of the stoichiometric matrix. 

Let us take the system given in the example of Eq.[T]{TJ and further decompose it into small reactons (atoms): 

CH 4 + CH 2 ^ C 2 H 6 2 (22) 

C 2 H 6 2 + 2 ^ H 2 2 + C 2 H 4 2 (23) 

C 2 H 4 2 ^ 2CH 2 0. (24) 



This system is only for the sake of demonstration and has no real counterpart. A realistic system will be considered in 
the penultimate section. X x = CH 4 0, X 2 = CH 2 0, X 3 = C 2 H 6 2 , X 4 = 2 , X 5 = H 2 2 , and X 6 = C 2 H 4 2 . 



There are three atomic reactons: A\ = C, Ai = H and A3 = O. The mass balance in each of these reactons can be 
seen by the fact that all reactions are equilibrated. The matrix writing of this is: 



a = [oi| 02! 03] 
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It can be easily verified that the relation *a? 



is satisfied. This simple decomposition of the molecules of the 



reaction network into their atoms provides a first possible solution of the left null space. 
Null Space Base and Optimal Reacton Decomposition 

This a priori decomposition only gives a collection of reactons, that may not describe the whole null space. It is 

possible to compute a base of the left null space of v by a Gauss-Jordan elimination (|9|), represented by the matrix 

< — ► 
a : 

NuTl (V) : V = . . . \a* s ] (26) 

all the al vectors that are found respect Eq.|2T]and thus result in one possible reacton decomposition. There is always 
a linear mapping from any set of reactons to the original set found by computing the null space. It is thus possible to 
derive any new and more convenient reacton decomposition from this preliminary set. 

Of course, several basis are possible. It is therefore important to find an "optimal" decomposition. As the er^ 
represent the number of reactons present in each molecule, large values of a^f. imply small reactons (because it means 
that more reactons are required to build the molecule). As a consequence, researching the linear combinations of the 
cjfe that compose the Null ( v ) while minimizing the individual and maximizing the number of nought values 
will lead to a desired maximization of the size of the reactons. Obtaining a base of the null space composed of such 
larger reactons present great interest, as the larger molecular subparts that are never broken through the reactions are 
likely to represent fundamental building blocks of the system. This should lead to an optimal network decomposition 
characterized by a small number of reactons. 

Computing the left null space base of the example system, by using the Gauss pivot method, gives the three 
following null vectors: 
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For a good representation of reactons, all the vector elements must be positive or zero integers. A new base can be 
searched for, by linearly combining the null vectors, aiming at maximizing the number of zero elements in each new 
vector. 

The second and third reacton possess only positive values and thus can be kept as such. In order to eliminate the 
negative values of the first one, we can just add the two first reactons: 
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This decomposition can be seen as a better one than the atomic decomposition. The molecules are at most de- 
composed into 3 reactons, instead of up to 10. Moreover, two molecules (X% and X4) are even composed of only one 
reacton, and can thus be identified as elementary building blocks of the system. 



Elementary decomposition of reactons 

Since we can express any reacton as a linear combination of the vectors of the null space base V, it is possible to 
link the elementary reactons to the atomic reactons. The elementary decomposition V = (oi| . . . \a£), for a system 
composed of a atoms from A\ to A a , can be written as: 
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The reacton decomposition V composed of s reactons from S\ to S s , can be written as: 
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Expressing a in terms of the null space base is equivalent to finding a matrix T so that: 
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Combining Eq. [29] and Eq. l30l gives: 
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Because a is full rank, it is possible to obtain the decomposition of the reactons into atom by: 

< > < X > 

S = T A. 

In the example of Eq. 12211241 linear combinations between al and can easily be found: 
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As a consequence, this new decomposition, aiming for bigger reactons, gives the three following ones: Si = O2, 
S 2 = H 2 and S 3 = CH 2 0. We thus have X 1 = (H 2 )(CH 2 0), X 2 = (CH 2 0), X 3 = (H 2 )(CH 2 0) 2 , X A = (0 2 ), 
X5 = (0 2 )(H 2 ), Xq = (CH 2 0) 2 . It can be checked that all the reactions of the example can be written as association 
and dissociation of these only three sub-elements. 



Transformation decomposition 



It is thus possible to find an automatic decomposition of molecules into reactons, and to find the atomic decomposition 
of these reactons into atoms, as long as the formula of molecules is known. 



The decomposition of molecules Xi into s reactons Sk can then be applied to the transformations themselves, 
allowing a focused study of each kind of reacton inside the system and how they do interfere. The molecules are 
decomposed as: 

s 

Xi = Y,<ri,kX$ k) - (40) 
k=i 

(k) 

At this point, it is important to track the position of reactons in the molecules they belong to. X- represents the 
subpart of a molecule Xi that corresponds to a reacton S^. 

We can decompose each reaction Rj into s partial reactions, each one describing the transfer of reactons from one 
compound to another: 
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Because there is a mass balance in each reacton Sk, we can decompose the reactions into partial reactions relative to 
each reacton: 
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Rj represents the subpart of the reaction Rj that involves only the reactons Sk. The symbol o represent the Hadamard 
product (i.e. elementwise product). 

The transformations relative to the reacton Sk can be summarized by the following stoichiometric matrix: 
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This operation correspond to multiplying each row by the corresponding element of the reacton vector. As each 
reacton vector is in the left null space of the stoichiometric matrix, we have: 



(47) 



=> T T V (fc) = 0". (48) 

That is the sum of the rows of each is null or, said differently, that the vector 1 is in the left null space of each 
one of these matrices. 

The system given in example can be decomposed into the following three reacton stoichiometric matrices: 
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A more compact notation can be used. Submatrices can be extracted for each reacton by removing the reactions 
that do not involve this reacton (i.e. columns) and the molecules that do not contain it (i.e. lines): 
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Given this simplified notation, we need to keep track of the meaning of each line and column for not losing the infor- 
mation about the involved reactions and compounds, but it definitely gives a simpler view of the reacton subsystem. 



Decomposition into fluxes and cycles 

Reactons are never broken into pieces, they are just transferred from one molecule to another. From this point on, all 
the reactions can be decomposed into simple A —>■ B transformation reactions. The whole system can be decomposed 
into either linear fluxes (a succession of transformation A B C — > • • • ) or cyclic fluxes (a succession of 
transformations A — ► B — > • • • — > A). This property of the new subsystems will greatly ease their analysis. 



Order of the reactions 



Because of the conservation of the reactons, the sum of the components of each colon is null (Eq. [48b - The sum of 
the positive numbers is thus identical in absolute value to the sum of negative numbers. This absolute value gives the 
number of reactons engaged in the reaction that is the order of the reaction for a given reacton. In order to linearize 
the system, it is necessary to determine the order of each reaction, so that each reaction of order n can be divided into 
n partial reactions of order 1. 

By defining the operations () + and ()~ as follows: 
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We can define the following function, giving the order of each reaction: 
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In our example, we obtain: 

Srd(V (1) ) = (0 1 0) ; Srd(V (2) ) = (l 1 0) ; ord (V (3) ) = (2 2 2) 
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Connection of the system to external fluxes 

Each combination *H of reactions leads to a global transformation t = 4r T^*c T . Provided this resultant transformation 
t can be compensated by an external flux / = — t , this combination can be maintained in a steady state. The 
stoichiometric matrix can then be extended by multiplying each reaction according to *c, and adding a new column 
representing this flux: 
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We now obtain matrices whose sum of components of each row (V 1 T and 1 T ) is null. 

Provided this addition of some external fluxes of reactons, the resulting system can be maintained in an active 
steady state. It is then possible to compute: 
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Each reacton turns out to be involved in a same number of creations and destructions. 



In the system given as example, we will allow exchanges of compounds X\, X2, X4 and X5, while compounds 
X3 and Xq will be considered as internal compounds. V is written as in Eq. [5] We first need to identify the vector 
c in order to guarantee the conservation of internal compounds. We also need to identify the vector t in order to 
guarantee the conservation of the external ones. Since X3 is only present in reactions R\ and R2, its conservation 
demands to combine a same number of R\ and R2. Xq is only present in reactions R2 and i?3, so that a same number 
of R2 and i?3 is required. As a consequence, there is only one possible combination of reactions that can lead to a 
steady state, V = (l 1 l). We then obtain: 
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And for each reacton: 
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Denoting by F the part of the matrix responsible for the external fluxes, the reduced expression of the matrix is: 
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Similarly, the reduced expression of this second reacton matrix is: 
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The reduced expression of the third reacton matrix is: 
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Cycle decomposition 

It is always possible to decompose the system into unitary transformations. All molecules can be decomposed into 
reactons and every reaction can be divided into partial sub-reactions concerned with the transfer of one reacton from 
one molecule to another. This amounts to an extension of the reacton matrix by duplicating each line and row according 
to their respective order. 

o 

For instance, the decomposition of V^ 3 ) (i.e. the matrix associated with the third reacton) is detailed in Fig.[2la 
and goes as follows. The elements of the matrix are distributed in such a way that only unitary reactions are obtained 
and every reacton is involved in only one reaction as a reactant and only one reaction as a product. The first step 
consists in supplying additional null rows and column according to the orders of molecules and reactions of the 
concerned reacton. In the second step, the matrix is parsed, first from left to right, then from the top to the bottom. 
If a non zero number is found in a column that already contains another number of the same sign, it is moved on the 
right. Then, if a non zero number is found in a row that already contains another number of the same sign, it is moved 
down. A number greater than 1 or smaller than — 1 must be considered as several 1 or — 1 occupying a same matrix 
cell. Consequently, this cell is adapted so as to contain only 1 and — 1 located as described above (no number of the 
same sign in a same row or column). 

A new matrix is thus obtained in which every line and every column contains one and only one pair of —1 and 
+1. As a result, the molecules have been decomposed into single reactons and the whole reaction network into single 

(k x) 

reacton transformations from one molecule to another. In this new matrix, the A - ' are associated with one of the 
reactons Sk of the molecule JQ. Since these apparently different reactons turn out to be identical, they can be swapped 



without altering the system. is the flux of reacton Sk, and R - are the partial reaction x of the reaction R - , 
corresponding to one single conversion of a reacton Sk from one molecule to another through the reaction Rj. Various 

(k x) 

combinations are actually possible, and the different R - ' can exchange their reactants or their products without 
altering the system. 

The right null space of this new matrix can simply be computed as the addition of single columns (Fig. |2]b). For 
a given reacton matrix, there is a simple way to proceed. Starting from a non zero value, one can go to the other 
non-zero value of the same line, then to the same non-zero value of the same column, then to the same non-zero value 
of the same line, etc. When the starting value is reached, the sum of all the visited columns is , and the sum of all 
the visited rows is . The following part of the right null space gives rise to a cycle: 

X™ J£l x<»-> ^1 Xt a) ^ X< 8 -> ^ Xt b) ^ Xt b) ^ X™ *^ (*<»>) (69) 

This one unique cycle involves all reactons and reactions. We can notice that, for example, the cycle goes twice 
through the molecule X3, to perform the reaction R2. Since the two reactons are actually identical inside this molecule, 
they can be exchanged. This amounts to the substitution: 

...Jh^x™ ^— + Xf' a) --- (70) 

R?' b) (3 b) Rf b) (3 b) 
1 — >Xf> — — > Xp b >--- (71) 

with: 

••• - ..Y : f" ; *2— > (72) 

...^_ x (3,6) X (3,a)... (73) 



The new version of the decomposition leads to two different combinations of reactions in the null space: 

[l | 1 | 1 I 1 ] = i?f o) + R { 2 ,b) + Rs ,a) + F (3) (74) 

[0 1 | 1 | 1 | ] = #f 6) + R { 2 ,a) + R { 3 ,b) (75) 

They are complementary (their sum is *c) and there is never more than one partial reaction involved. The system is 
thus completely decoupled, keeping apart all particular reactive fluxes. If there were still crossings like the one just 
resolved, a new matrix should be formed by similarly inverting the crossing point. This process can be continued 
as long as crossings still exist, until a fully decomposed reaction network is reached. With such a decomposition 
procedure, the following cycles are obtained: 

X (3) *tl X m #t X™ ^1 X™ ^ (xf ) (76) 

x (3,a) x (3,a) «^ ^3,a) «^ (x (3,a) } (77) 

Moreover, actually represents an incoming flux of xf (noted -Ff ) and an outgoing flux of xf (noted 
F_ (3)), so that the cycle of Eq. l76l represents a linear flux, which can be written: 

7?( 3 ' , v d(3><0 , v of 3 ' 6 ) , D (3,f)) , , P (3) 

xf xi a » X<M X™ ^ (78) 

These two fluxes are coupled since the involved partial reactions are always working together (e.g. i?f = i?f + 
R { *' b) ), and the two pairs of reactons {xf > a) , xf ' b) } and {xf ' a) , xf ' 6) } are linked into the same molecules X3 and 
Xg. As X2 is only composed of one S3 reacton, xf a ^ and xf ^ are not linked, but represent two reactons that are 
present in two different X2 molecules. This can be represented in a graphic form (Fig. [3^) or, more simply, by just 
emphasising the fluxes (Fig. EJ?). The linear flux and the cycle are coupled by the three reactions. 

The treatment of the two other reactons is much simpler. In each case, only one linear flux is observed. 

For the first reacton: 

-^U xf ^-^ xf ► (79) 



For the second reacton: 



^Af'^xf^xf^ (80, 



Rebuilding the complete system 

According to ord (^^j (Eq.[59l), the system is composed of lXi + 2X2 + IX3 + IX4 + IX5 + IXq. These molecules 
are decomposed as follows into the three reactons: 





= xf +xf 


(81) 


xf 


= xf ^ 


(82) 


xf 


= xf fe) 


(83) 


x 3 


= xf + xf a) + xf b) 


(84) 


x 4 


= xf 


(85) 


x 5 


= x® + xP 


(86) 


x 6 


= xf a) + xf ,b) 


(87) 



With respect to the flux analysis, the system is represented by V = (l 1 l), that is IR\ + IR2 + I-R3. These 
reactions are decomposed into the following partial reactions: 

Ri = 4 2) +4 3 ' 0) +4 3,6) (88) 

R 2 = i?f+4 2) +4 3 ' a) + 4 3,6) (89) 

Rs = i?f (90) 



The whole system can then be rebuilt by associating the corresponding reactons and the partial reactions, leading 
to Fig. [3fc. The flux of Si is represented in red, the flux of S2 in green, and both the flux and cycle of S3 in blue. 
The dots indicate the reactions where the fluxes are coupled. The link between the fluxes represent the molecules, 
composed of several reactons. 

The total system has thus been decomposed in a way that emphasizes the evolution of its different subparts. There 
is a global flux from X\ to X2, another one from X4 to X5, and an internal cycle of reacton S3. The autocatalytic 
property of this system can be seen by the coupling of a linear flux and an internal cycle concerning the same reacton 
S 3 . 

Application to a Realistic Network 

In this section, a more complex and realistic example, a partial E.Coli metabolism (10) is treated by applying the 
same sequence of algorithmic operations. This system describes the decomposition of glucose into carbon dioxide. 
It involves the transformation of 37 molecules through 28 reactions. We have kept the same notations as described 
by Beard et al. (llOh . except for HEXT (exchange of H + through a membrane) that has been replaced by two 
compounds HI and H2, corresponding to internal and external protons, in order to keep the mass balance. This 
system is described in Fig. [4] and the corresponding stoichiometric matrix in Fig. [5] 

A base of the left null space can be computed from this matrix, leading to the molecular decomposition given in 
Fig. The molecules are reduced to the combination of the following reactons: 



Si 


= NADH 


(91) 


S-2 


= COA 


(92) 




= Pi 


(93) 


Si 


= NADPH 


(94) 


s 5 


= FADH 


(95) 


s 6 


= AMP 


(96) 


s 7 


= O 


(97) 


S 8 


= QH 2 


(98) 


So 


= SUCC 


(99) 


S10 


= C 


(100) 


Sn 


= H 


(101) 



We must note that, in the equations, the water molecules are implicit, and thus do not appear in this decomposition. 
This explains why the glucose (X3) is written as 6S10 = C% rather than C6H12O6 = Ce(H20)6. 

We can see here how the 37 chemical compounds of the network can be reduced to a combination of only 11 
reactons. If this decomposition is an obvious one for any biochemist, it is important to understand here that it can 
be automatically obtained, with no further knowledge that the stoichiometric matrix. On the basis of this reacton 
decomposition, it is then be possible to focus on the evolution of some given reacton, e.g. the evolution of C in 
the metabolism, from glucose to carbon dioxide, the evolution of O from dioxygen to carbon dioxyde, the use of Pi 
throughout the whole network, etc. 

For example, following the reacton £10, that is the carbon coming from glucose, we can reduce the whole stoi- 
chiometric matrix V* of dimension 37 x 28 to the stoichiometric matrix of the sub-network relative to S10 V^ 10 ) of 
dimension 16 x 16 (see Fig. IT]). It can then easily be decomposed into single fluxes and cycles (see Fig. [8]). It becomes 
more tractable to identify the progressive degradation of glucose, each carbon following a linear flux towards carbon 
dioxide and being released in three possible places. The whole system is coupled to the PEP/PYR cycles. 

Conclusion 

We have shown in this article that the left null space analysis of the stoichiometric matrix can lead to an automatic 
decomposition of molecules into physically meaningful sub-elements called "reactons". Besides giving insight to the 
different moieties that can be studied through the network, the discovery of reactons leads to a natural simplification 



of the network, by dividing it into subnetworks, each one related with one specific reacton. These subnetworks can 
be easily studied and understood in terms of simple fluxes and loops, by separating the reactions into unitary partial 
reactions, describing the transfer of one reacton from one molecule to another. The global network can then be seen 
as a coupling of these elementary sub-elements. All these algorithmic manipulations can be grouped into one single 
software that remains simple to implement and use. Once the reactons have been identified and the corresponding 
decomposition into subnetworks achieved, the different modes are readily obtained. 

However, the simplification allowed by such a decomposition into reactons is nevertheless offset by the difficulty 
of deriving an optimal reacton decomposition. This amounts to computing a sparse null space base, which is far from 
being a trivial problem (EI 03). This step remains however of fundamental importance, as the sparsest null space will 
lead to the largest interesting reactons and the simplest corresponding sub-networks. The "brute-force" approach - i.e. 
computing a null space from the classical Gauss-Jordan elimination (9) - is easy to implement, but only leads to an 
optimal solution following a huge amount of possible linear combinations of vectors. This approach is not realistic for 
large systems on account of the exponentially increasing number of operations. This problem has nevertheless been 
thoroughly studied in the literature (11 lL Il2l Il3l Il4) . and a careful examination of such work should help the future 
development of new algorithms that are better adapted for obtaining the optimal reactons. 
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Figure Legends 



Figure [TJ 

Representation of the example chemical network given in Eq. [JJ|3] 
Figured 

o 

Decomposition of *z7^ 3 ). a) Detail of the operations leading to a square matrix representing unitary transfers of 
reactons from one molecule to another, b) Detail of the operations leading to the decoupling of cycles into elementary 
independent cycles. 

Figure |3] 

Graphical representation of the flux/cycle decomposition of the example chemical network of Fig. [TJ a) Complete 
decomposition of the flux and cycle relative to the reacton S3. The partial reactions are represented in red, and the 
reactons are represented in violet, b) Simplified representation of the flux and cycle relative to the reacton £3. The 
dots represent the coupling between partial reactions inside the complete reactions, c) Flux/cycle decomposition 
representation for the whole network. The flux of Si is in red, the flux of S2 in green and the flux and cycles of S3 in 
blue. The dots represent the link between partial reactions, and the segments represent the link between reactons. 

Figure S] 

Chemical network representing a partial metabolism of E.Coli (10), composed of 37 molecules and 28 reactions. 
Figure H 

Stoichiometric matrix of the E.Coli chemical network described in Fig. [4] 
Figure |6] 

Sparse base of the left null space of the stoichiomeuic matrix of the E.Coli chemical network of Fig. [5] 
Figure |7] 

Matrix representations of the subnetworks of the E.Coli chemical network, relative to the flow of organic carbon 
(reacton Sio). 

Figure! 

Flux/cycle decomposition of the subnetworks of the E.Coli chemical network, relative to the flow of organic carbon 
(reacton Sio)- The dots represent the link between partial reactions, and the segments represent the link between 
reactons. *: Reactions R\, R2 and R3. f: Reactions 2Rt, 2R$ and 2R$. p. Reactions 2Ri§ and 2R\q. 
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Figure 5: 
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